Sweeping a molecular Bose-Einstein condensate across a Feshbach resonance. 



o 
o 

(N 
-i— > 
o 

O 

in 

(N 



5-1 



o 
o 



> ■ 

00 ' 

-xf : 
\o . 
o ■ 

: 

o . 
-i— > ■ 

I ' 
T3 ■ 

c : 
o . 

o ■ 



X 



Masudul Haque* and H. T. C. Stooff 
Institute for Theoretical Physics, Utrecht University, 
Leuvenlaan 4, 3584 CE Utrecht, The Netherlands 
(Dated: February 2, 2008) 

We consider the dissociation of a molecular Bose-Einstein condensate during a magnetic-field 
sweep through a Feshbach resonance that starts on the molecular side of the resonance and ends on 
the atomic side. In particular, we determine the energy distribution of the atoms produced after 
the sweep. We find that the shape of the energy distribution strongly depends on the rate of the 
magnetic-field sweep, in a manner that is in good agreement with recent experiments. 



Introduction. — Feshbach resonances have established 
themselves as one of the most important tools in the 
area of ultracold atomic gases [1]. These resonances pro- 
vide the opportunity to make use of an applied magnetic 
field to tune the interactions between the atoms, essen- 
tially by changing the energy of a diatomic molecular 
state from a region where the molecules are more sta- 
ble to a region where the atoms are more stable. For 
bosonic atoms, Feshbach resonances allow for coherent 
oscillations between atomic and molecular Bose-Einstein 
condensates [2], Moreover, in a Fermi gas these Fesh- 
bach resonances recently have made a novel BCS-BEC 
crossover phenomenon accessible to experiments [3-7] . 

In this Letter, we study the response of a spatially 
uniform molecular condensate to a magnetic-field sweep 
across a Feshbach resonance. The essential physics of the 
nonequilibrium dissociation of the molecules is the same 
for bosonic and fermionic atoms, except for the feedback 
effects arising from the atomic statistics. We show that 
under appropriate conditions these feedback effects can 
be neglected, and in that case our results are equally valid 
for fermionic and bosonic atoms. 

We focus in particular on two quantities. First, we de- 
velop a method for determining the decay of the molec- 
ular condensate during the sweep. Second, we investi- 
gate for different kinds of sweeps the resulting energy 
spectrum of the dissociated atoms, which turns out to 
be closely related to the time evolution of the molecular 
condensate. Our approach to the time evolution is based 
on the equation of motion for the molecular condensate 
wave function <j) m (t), developed by Duine and Stoof [I, 8]. 
While no experiment has yet probed the molecular disso- 
ciation as a function of time, the decay of the molecular 
condensate density is in principle experimentally accessi- 
ble by a series of destructive measurements. 

Measurements of atomic energies after a dissociation 
ramp have been performed recently by Mukaiyama et al. 
[9] and by Diirr et al. [10]. The ramps of Ref. [9] are 
restricted to cases where all the molecules dissociate dur- 
ing the sweep. The shape of the dissociation spectrum 
for this type of sweep is well described by an analyti- 
cal expression, derived by Mukaiyama et al. [9] and also 
by Goral et al. [11]. For such a relatively slow type of 
ramp, our calculated spectra agree with this analytical 



form as well. In Ref. [10] also a different, more rapid, 
type of sweep experiment is reported, in which most of 
the dissociation occurs after the sweep is over. As a 
result there is now a sharp peak in the atomic energy 
distribution corresponding to half the molecular energy 
at the end of the sweep. These spectra are not described 
by the above mentioned analytical result. Our approach, 
however, allows us to calculate also the spectra for such 
fast sweeps. 

Both these sweep experiments involve bosonic atoms, 
in particular 23 Na and 87 Rb, respectively. An additional 
motivation for investigating sweeps in the molecule-to- 
atom direction comes from considerations related to the 
cooling of fermionic atoms such as 6 Li and 40 K. To avoid 
the problems due to Fermi statistics, one proposal is to 
first cool diatomic molecules to Bose-Einstein conden- 
sation, and then to slowly sweep the system across the 
resonance to the atomic side. However, the effect of the 
sweep rate on the resulting temperature of the atomic gas 
is not well-understood at present. It is, therefore, impor- 
tant to study the dissociation process and the spectrum 
during a sweep across a Feshbach resonance in some de- 
tail. 

Generalized Gross-Pitaevskii equation. — A Feshbach 
resonance is characterized by its magnetic field location 
B , its width AB, the background scattering length a^g, 
and the magnetic-moment difference A/i between two 
atoms and the bare molecule associated with the Fes- 
hbach resonance. Assuming that these quantities are for 
the atomic gas of interest known from experiment, the 
equation of motion for the bare molecular condensate 
wave function <p m (t) becomes [1] 



ihd t 4> m (t) = S(t) - ij]\Jihd t 4> m (t) 



(1) 



where 5(t) = A/j,(B(i) — Bo) is the so-called detuning 
that is time dependent during the sweep. The imagi- 
nary term in this equation of motion describes the decay 
of molecules into pairs of atoms and is determined in 
the Bose case by the quantity 77 = m 1 / 2 ab g Ai?A/z//i = 
m 3//2 g 2 /27rfi 3 , where g is the atom-molecule coupling. 
When the atoms are fermionic, the right-hand side is a 
factor of two smaller. Note that rj 2 has the units of en- 
ergy and physically quantifies the width of the resonance. 
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The imaginary term arises from the molecular self energy 
HZ m {z)= =^<^>^= =^V~z, (2) 



where the single and double lines represent atomic and 
molecular propagators, respectively. The square root 
of energy reflects the three-dimensional final density of 
states expected from a Fermi's Golden Rule calculation 
of the dissociation process. The \/ihdt form in Eq. (1) 
arises from evaluating the complex self energy in Eq. (2) 
from above on the real axis. 

The above self energy leads physically to a dressing of 
the molecules and as a result to a molecular density of 
states that for negative detuning contains a delta func- 
tion with weight Z = (1 + r)/2y/\e m \)~ 1 at the dressed 
molecular binding energy e m (<5) = 8 — rj 2 /2 + r]y / 'ri 2 /A — 8 
[1]. Since the wave function (f> m appearing in Eq. (1) rep- 
resents the condensate wave function for bare molecules, 
the dressed molecular density is a factor Z~ x larger than 
the bare molecular density and equal to \<fi m \ 2 /Z. More- 
over, the bare molecular density can now change in two 
ways. Coherent changes involve a shift in the linear su- 
perposition making up the dressed molecular wave func- 
tion, i.e., a change of the bare molecular fraction Z ', while 
the dressed molecular density remains constant [1, 2]. In- 
coherent changes are associated with the break-up of the 
dressed molecule itself. 

The generalized Gross-Pitaevskii equation in Eq. (1) 
involves some approximations. First, we are considering 
sufficiently low molecular densities that the molecule- 
molecule interactions can be neglected for most of the 
sweep duration, which is justified for the experiments of 
interest to us. Second, we are not including the effect 
of the thermal part of the molecular gas. This is justi- 
fied for sufficiently low temperatures. Most importantly, 
we are neglecting feedback effects, i.e., the dynamics and 
statistics of the atoms that are formed during and after 
the sweep. In particular, we do not allow for their re- 
association. The neglect of atomic statistics, i.e., Bose 
enhancement or Pauli blocking, is reasonable for mono- 
tonic sweeps, because the dissociation then always occurs 
at different energies. In that case we can estimate that 
the atomic statistics can be absolutely neglected if the 
molecular density obeys n m < M / \f2it 2 g 2 . Not surpris- 
ingly, this approximation is better for faster ramps, and 
worst for a constant positive detuning. 

We now outline the method of solving Eq. (1), given 
a particular time-dependent detuning 8(t). We are espe- 
cially interested in magnetic-field sweeps across the reso- 
nance that take the detuning 8 from some negative value 
8\ at t = to the final positive value <5f at time t = T. 
Introducing an auxiliary constant detuning 8 C , we can 
write 



The solution of this differential equation is the sum of a 
particular solution and the solution of the corresponding 
homogeneous equation. We thus have 



df G(t ~ t') [s(t') ~ <y </> m (0 



(4) 

Here G(t) is the Green's function for the operator acting 
on 4> m (t) in the left-hand side of Eq. (3), i.e., 



G(t) 



2tt tujj+ — rj\/—fku+ — S c 



(5) 



where fku + = hut + iO. The real part of the pole in G(oj) 
gives the molecular energy and the imaginary part gives 
the dissociation rate, both at the detuning 8 C . More- 
over, ip(t) is the solution of the homogeneous equation of 
Eq. (3). For 8 C < 0, the solution is an oscillatory func- 
tion, ij){t) oc e~ JWct / ri , with a frequency determined by the 
molecular binding energy corresponding to the detuning 
(5 C , i.e., hw c = e m (S c ). In the following we use 8 C = Si, 
because the integral in Eq. (4) then has a lower limit of 
instead of — oo. 

The retarded Green's function G(t) is calculated by 
closing the contour in the negative half-plane. The inte- 
grand has a branch cut because of the square root in the 
energy denominator. For S c < 0, there is also a pole at 
e m . Taking these into account we obtain 



G(t) = 9(t) 



IT] 
7T 



dtu- 



huje 



-iuit 



(huj — S c ) 2 + rj 2 haj 



rj + 2\/—Huj c 



ihdt 



\fihd t 



</» m (i) = [S(t)-S c ]<p m (t) ■ (3) 



Having calculated the Green's function G(t) and the wave 
function ip(t) corresponding to S Cl we can numerically 
solve Eq. (4) for (f> m (t). 

Magnetic field sweeps. — In Fig. 1 we show the de- 
cay of the bare molecular wave function <fi m (t) for a par- 
ticular ramp that is also shown in the top panel. The 
center panel shows the oscillatory decay of the real part 
of 4> m (t). The imaginary part yields similar qualitative 
information and is not shown. The dotted background is 
the nondecaying oscillatory 4>m{t) for a negative detun- 
ing that is held fixed at 5i. We see that <j) m {t) for the 
ramp has a gradually decreasing oscillation frequency, 
corresponding to the decreasing magnitude of the bind- 
ing energy e m (<5) as 8 approaches zero from the negative 
side. The bottom panel shows the decay of the bare den- 
sity n m (t) = \<j) m (t)\ 2 . The decay for S(t) < is mostly 
coherent, so that the dressed molecule density remains 
almost fixed. On the other hand, the decay at positive 
detuning is incoherent, and the decay of the bare density 
n m (t) corresponds to the decay of dressed molecules. 

The ramp in the experiment of Mukaiyama et al. [9] 
involves an initial and final detuning £j.f of about =p50ry 2 , 
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FIG. 1: Decay of the molecular wave function m (i)/|(/> m (O)| 
with time, when the detuning changes linearly with time from 
Si = — IO77 2 to 5f = +IO77 2 as shown in the top panel. The 
center panel shows the real part of (j> m (t). The bottom panel 
shows the decay of the molecular density, n m (i)/n m (0) = 

|<Mt)| 2 /l<Mo)| 2 . 



and ramp times of T > lOOh/r/ 2 . For the results shown 
in Fig. 1, we used smaller values for both, so as to have 
fewer oscillations to visualize. 

We now outline how the incoherent dissociation spec- 
trum is obtained, once <fi m (t) has been calculated. We cal- 
culate the time derivative of the molecular density from 
the equation of motion in Eq. (1). Not surprisingly, only 
positive frequencies contribute, because the molecules are 
stable for negative energies. We obtain 



dt 



incoh 



v_ 

irh 



dwy frhj |(/) m ([j)| , (6) 



where 4> m (tjj) is the Fourier transform of </> m (£). The inte- 
grand in the right-hand side is interpreted as the molec- 
ular dissociation spectrum f mo \(e) = r\\[e |</> m (e)| 2 /irh. 
Since each atom has half the energy of the dissociating 
molecule, / mo i(e) corresponds to the atomic spectrum at 
e/2, i.e., / a t(e) = 4/ mo i(2e). The additional factor of 2 is 
because each molecule gives rise to two atoms. 

Before presenting results for the spectrum, we briefly 
discuss the analytical result mentioned previously [9, 11]. 
This result can be derived by assuming that for a par- 
ticular detuning, each decaying molecule has a sharply 
defined energy. This means that the width of the molec- 
ular spectral function is neglected. Moreover, it is implic- 
itly assumed that all the molecules are dissociated before 
the sweep is over. If we denote by n m (e) the density of 
molecules at energy e and use dt = de/e ~ de/5, we have 




e/r| 

FIG. 2: Spectra for cases where the analytic derivation is 
valid, because the molecules decay almost completely before 
the sweep is over. Insets show the detuning ramps with differ- 
ent ramp rates, and the bare molecular condensate densities. 



n m (e + de) — n m (e) — — (27]y/e/H5)n m (e)de. Solving for 
n m (e), the molecular dissociation spectrum follows then 
from / mo i = dn m (e) / de. The atomic energy spectrum / at 
can be deduced from / mo i as before with the result 



/at(e) 



n m(0) — —7 — exp 
no 



47? (2e) 



3/2 



3M 



(7) 



We stress that this expression is not a Landau-Zener or 
adiabatic result. Its derivation does not require the sweep 
to be slow, although it is required that the sweep con- 
tinues until all the molecules have dissociated. Moreover, 
the derivation relies on the approximation e ~ <5, which is 
not accurate at small energies e <C i] 2 , for which we have 
that e ~ 5 2 /rj 2 . In Fig. 2 we show calculated spectra for 
two cases where the above mentioned conditions holds. 
The spectra indeed have the exponential form of Eq. (7), 
however they are in general normalized according to the 
initial dressed molecular density n m (0)/Z; instead of the 
initial bare molecular density n m (0). The ramp repre- 
sented by the solid curves has a higher ramp rate and 
hence the resulting spectrum is tilted to higher energies. 
Note that experimentally, only the total atomic energy 
has been measured as a function of ramp speeds [9, 10]. 
These measured total energies are fit well by the total en- 
ergies calculated analytically as J de e / a t (e) . Since our 
spectra match the analytic result for the experimentally 
relevant case of Zi ~ 1, agreement with the measured 
total energies follows, and is not shown here. 

In Fig. 3 we show spectra calculated for two cases 
where Eq. (7) does not hold because the sweep is over 
before all the molecules have decayed. This type of sweep 
has been employed in a recent experiment [10]. The den- 
sity images in that experiment indicate that the spectrum 
is peaked at the energy corresponding to the final detun- 
ing. This feature is clear in Fig. 3. The dashed curves 
represent an infinitely fast ramp. The spectrum shows a 



FIG. 3: Spectra for cases where the analytic form is not valid. 
Dashed curves: Infinitely fast ramp at which the detuning 
switched from Si = —6r/ 2 to 5{ = +6w 2 at t — 0. Full curves: 
Fast ramp that is over before all the molecules can decay. 

peak around 5{/2, but is spread over a relatively large 
range of energies. The spread demonstrates that our 
method takes into account the distribution of molecular 
energies, i.e., the width of the molecular density of states 
at a fixed positive detuning. The solid curves in Fig. 3 
represents a ramp that is over before all the molecules 
have decayed. A significant fraction decays while the de- 
tuning is constant at its final value. The atomic spectrum 
here also shows a peak at an energy around half the final 
detuning. 

In order to demonstrate mainly coherent loss of bare 
molecules, we display in Fig. 4 a sweep that ends at neg- 
ative detuning, i.e., before reaching the resonance. The 
density of bare molecules decreases, but most of this de- 
crease is due to the coherent change of the dressed molec- 
ular wave function. The final density m is slightly smaller 
than (Z{ / Zi)n m (0) , where Zi and Z\ are the weight of the 
delta function in the molecular density of states at the 
initial and the final detuning, respectively. The difference 
between rif and (Z{/Zi)n m (0) appears as the total weight 
of the incoherent spectrum. Note also that the density 
undergoes some coherent oscillations before settling at 
its final value. These are internal Josephson oscillations 
of the dressed molecular condensate. 

Conclusions. — We have shown how to calculate the 
time evolution of the molecular condensate wave function 
cj) m (t) during sweeps across a Feshbach resonance. We 
have also shown how the time dependence of this wave 
function leads to the dissociation spectrum. Our calcu- 
lation for the atomic energy spectrum is more general 
than the result obtained by Mukaiyama et al. [9] and by 
Goral et al. [11]. The theory presented here is applicable 
for ramps of essentially arbitrary shape and initial condi- 
tions. This has been demonstrated with the fast ramps, 
in Fig. 3. Moreover, we have obtained the, in first in- 
stance surprising, result that a nonzero dissociation spec- 
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FIG. 4: A sweep that starts and ends at negative detuning. In 
this case, the final density is expected to be (Zf/i?i)n m (0) ~ 
0.7775n m (0). The actual final density, nt ~ 0.735n m (0), is 
smaller, in agreement with the nonzero atomic spectrum. 



trum can also be obtained by a changing magnetic field 
that always remains below the Feshbach resonance. 

Finally, we point out two important issues that remain 
unsolved. First is the problem of incorporating the effect 
of the atomic statistics. In principle, this can be achieved 
by coupling the generalized Gross-Pitaevskii equation for 
the molecular condensate to a quantum Boltzmann equa- 
tion for the atoms. It remains to be seen, however, if this 
leads to a feasible approach. Second, our calculation is 
based on the approximation for the molecular self en- 
ergy given in Eq. (2). For broad Feshbach resonances, 
such as the one in 6 Li that is widely used at present for 
the study of the BCS-BEC crossover [4-7] , the molecular 
self energy can no longer be expressed in terms of the 
r\ parameter alone, and depends also on the background 
scattering length abg [1] • Incorporating this improvement 
is also left for future work. 
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